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Abstract 

We study the influence of some quenched disorder in the sequence of 
monomers on the entropic elasticity of long polymeric chains. Starting from 
the Kratky-Porod model, we show numerically that some randomness in the 
favoured angles between successive segments induces a change in the elonga- 
tion versus force characteristics, and this change can be well described by a 
simple renormalisation of the elastic constant. The effective coupling constant 
is computed by an analytic study of the low force regime. LPTENS 97/51. 
New tools have recently been developped for the manipulation of single molecules, in 
particular DNA |R|3|,]3]J2] , protein fibers, such as titin |5|-|7| and polymers ||. In these ex- 
periments, the fibers are stretched by various means (optical or magnetic tweezers, flexible 
microscopic cantilevers, stokes drag) and their extension is measured. The various models 
(e.g., the Freely Jointed Chain (FJC) or Worm-like Chain (WLC) models) used to analyse 
the data are based on the elastic theory of an homogenous polymer stretched in its entropic 
regime. Because the actual molecules, particularly proteins and DNA, are in fact heteropoly- 
mers, it is interesting to study how their intrinsic quenched randomness modify their elastic 
behaviour ||. This could be particularly relevant to possible stretching experiments on 
single strand DNA or proteins in denatured conditions. 

In this paper we wish to understand the qualitative new features introduced by random- 
ness. We shall not consider a particularly realistic type of disorder (this should be tailored 



for each type of molecule). Rather we want to start from the usual Kratky-Porod model [10 
of an elastic chain, and compute the effects of the introduction of a certain simple type of 
disorder, namely some preferred orientations between successive links. In the following we 
compute analytically the effective (disorder averaged) persistence length at low extensions 
and use numerical methods (transfer matrix and Monte-Carlo) to solve for the behaviour 
of the chain at all extensions. Surprisingly, it appears that the full elastic response of this 
random chain is similar to the response of an homogenous chain with a different elastic con- 
stant. The effective elastic constant can be determined from the calculated low extension 
effective persistence length. 
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In the Kratky-Porod model, the energy 8 of an homogenous polymer chain, consisting of 
iV segments of size b (maximal extension L = Nb) and orientation given by the unit vectors 
ti is equal to: 
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where oii is the angle between the successive segments i — 1 and i and Q{ is the angle 
between segment i and the z-axis, the direction of the stretching force F. One recognises 
the analogy between polymer chains and a magnetic system: the ferromagnetic chain of 
vector spins in a magnetic field [EI]. The magnetic field H is related to the applied force 



F through H = Fb/ksT, while the ferromagnetic coupling K is related to the stiffness or 
persistence length of the chain £t by K = £r/b. The FJC model (obtained by setting K — 0) 
is equivalent to the behaviour of an ensemble of non-interacting spins. The WLC model is 
the continuum limit of Eq.(|l]): b — ► 0, £t and F finite (i.e.: K ^> 1 ^> H): 
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As the only non-dimensional parameter in Eq.(Q) is F^r/fc^T, the relative extension 
x =< cos 9 >= l/L (equal to the magnetization density in the spin language) will be a 
function of that parameter only, or inverting: 

F = (k B T/£ T ) g(x) . (3) 

The function g(x) (valid for any homopolymer) can be computed from a Hamiltonian formu- 
lation by solving for the ground state of a quantum dumbel |12| . It is known to reproduce 



with good accuracy the experimental data on DNA [T3| . 

In the following we shall consider a random version of the Kratky-Porod Model just 
described, where the energy of the chain is: 



k B T 



Here ipi is a preferred random orientation between the successive segments i — 1 and 
i along the chain. The angles ipi are quenched, independent and identically distributed, 
random variables, with probability distribution P(ipi). It is easy to obtain the spin-spin 
correlation function at zero field UTTl: 



i+l 

< ti.t i+l >=< ti.u+i^ > f{^i + i) = n f(ipj) , (5) 

j=i+i 

where: 

J da sin a cos a e Kco < a -^) 
/W = J da sin a e Kc °< a ~^ ' ^ ' 

The correlation (Q) decays at large argument as exp(— rb/£ e ff), where the effective cor- 
relation length is 6/^ e // = — log(/) (throughout the paper, the overline denotes the average 
over the distribution P{ip) of the random angle ip). This result for the effective correlation 
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length can be simplified in the continuum limit where 6^0. In this limit the stiffness scales 
as K = fr/b- We shall make the reasonable assumption that the distribution of disorder 
scales as P(ip) = P(ifj)/Vb where if> = ip/Vb, and that the angles ifti are all positive. Using 
the fact that the angle a also scales as Vb in @, one gets to first order in b: 

fW = 1 - ^G(y) (7) 

where y = y/£rif> an d the function G is: 

G(y) = 3 + y 2 - (l + ^| ye^ 2 (l + Erf [y/^2])) * (8) 
In terms of this function G the effective persistence length is given by: 

T- = VT [ d ^ H4>)G$>/b) (9) 

In the case where the disorder is small, which means that the typical values of the angles 



ipi are much smaller than y 6/£t, one can get a general result for the effective persistence 
length by expanding the function G at small arguments: 

7^ + \/??+f 1 -T^+- ( 10 ) 
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In the reverse case of strong disorder one gets: 
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In fig. 1 we plot the effective persistence length in units of the persistence length of the 
pure system, £ e ff/fr, as a function of the ratio r = £r/£d {t is a measure of the strength of 



disorder), for a half-gaussian distribution of disorder P(if>) = y 2^ d / Tib exp(—^ d ip 2 / 2b) (if) > 
0). The previous asymptotic results show that Ce/z/Cr behaves as 1 — a/t/2 at small r and as 
2/r at large r. These results differ from the heuristic arguments often found in the litterature 
nSyiEfl which appear to consider the annealed disorder case with gaussian disorder, yielding: 
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or Ceff/Cr — 1/(1 + t) , which is also shown in fig. 1. 

Very often, the actual measurements do not have a direct access to the correlation length, 
but they deduce it from low force measurements. In the limit of small forces, the relative 
extension is linear in the force with a slope: 

In principle, this quantity depends on the average of f(if)), and is not directly related to 
£eff (which depends on log(/)). This property of non self-averageness of the correlation 



function has been well studied in the litterature on disordered spin chains [TJJ. However it 
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Effective persistence length induced by the disorder 
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FIG. 1. The ratio of the effective persistence length to the one of the pure homopolymer, 
£eff/£,T, versus the strength of the disorder, measured by the ratio r = ^r/^d, in logarithmic scale. 
The full curve is the exact result from (||) for a half gaussian disorder (see text) , the dashed line is 
the result of the simple heuristic formula (|l2|). 

turns out that in the continuum limit b — > this subtlety disappears: using (0), we find 
that the slope is related to the effective persistence length through the usual WLC formula 
S = 2£ eff /(3k B T). 

The analytical results just presented are limited to the low force (linear response) regime. 
To study the elastic behaviour of a stretched heteropolymer for large extensions (x ~ 1), we 
will compute the partition function and the average extension (magnetisation) by transfer 
matrix and Monte-Carlo (MC) methods. 

The transfer matrix T(tj,tj_i) can be read from Eq.([|): 

T(ti,ti_i) = exp[K cos(ai — + H (cos 9i + cos 9i_i)/2] , (14) 

with: 

cos(cti - ipi) = cosipi ti ■ ti-! + smipi \Jl-(ti- fi_i) 2 

£j ■ £j_i = cos 6i cos 6i_i + sin 0; sin cos(0j — fa-i) ■ (15) 

The probability density p% of the angles 9i,(f>i obeys the recursion relation: 

Pi(0i,&;0o,0o) = C J waOi-xdOi^dtpi-i T(t h t^i) Pi-i(0i_i,&-i;0o,0o) , (16) 

where C is a normalisation constant. The relative extension (magnetisation) of long chains 
is given by x = liniAr^oo Tr[cos pj^}. By discretising the angular variables {9 i} (pi} and 
iterating Eq. (W&j , it is then straightforward to calculate numerically the relative extension x 
of a random chain due to a force F. In the limit of long chains this characteristics becomes 
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Random Kratky— Porod model 
f T /f d =0.1 ; £ T = 50 nm 




relative extension 

FIG. 2. The force elongation characteristics of the random chain, obtained by transfer matrix 
computations. The full lines are the predictions of the model without disorder, with its original 
elastic constant (lower curve), and with an effective elastic constant given by ® (upper curve). 

independent of the particular realization of the disorder. The results are presented in Figs. 
2 and 3, for the case of half gaussian disorder with two different values of r. The value of 
b/^T = -05 has been chosen small enough that the results are close to those of the continuum 
limit (in the absence of disorder the relative difference is less than 8 %). The continuous line 
is the force versus extension curve of a homogenous WLC, Eq.(|3]) with persistence length £ e ff 
given by Eq.(^). The remarkable and surprising agreement between the numerical data and 
the WLC results implies that for all experimental purposes the elongation curve of a random 
Kratky-Porod chain is identical to that of a homogenous chain with a smaller persistence 
length, i.e. it follows Eq.(|3]) with £ e // replacing £t- 

To further check this surprising result, we have run a Monte-Carlo simulation of a dis- 
cretised random Kratky-Porod chain with 6/£t — -05; N = 600 ; {ksT = 1). Rapid 
equilibration is achieved by a MC move which involves a possible global rotation of the 
chain for all segments j > i with i chosen at random). If the energy of the chain after the 
move is lower than before it [E new < £ id) the move is accepted with probability p = 1, 
otherwise it is performed with probability p = exp[— (£ new — £ id)]- The total number of 
Monte Carlo steps was around 2. 10 6 .In the absence of disorder (l/£d = 0) onr recovers the 
results of an homogenous WLC, Eq.(|3|). In the presence of disorder (£r/£d = 2) and after 
averaging over Md = 80 configurations of disorder, one recovers the transfer matrix results 
for the extension of a random Kratky-Porod chain (see Fig. 3). 

In summary, we have shown that a random Kratky-Porod chain with a certain type of 
disorder is well aproximated by a pure chain with an effective elastic constant. It would be 
interesting to check further the precision of this result, and its dependance on b, to see if 
it holds for wider classes of disorder, and to understand better why there is such a form of 
universality. 

We wish to thank R. Zeitak for useful discussions. D.B. acknowledges the partial support 
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Random Kratky— Porod model . 
f T /f d = 2.0 ; £ T = 50 nm 
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FIG. 3. The force elongation characteristics of the random chain, obtained by transfer matrix 
computations, and by Monte Carlo simulations. The full lines are the predictions of the model 
without disorder, with its original elastic constant (lower curve), and with an effective elastic 
constant £ e // given by ([!]) (upper curve). Two types of disorder, with the same £ e //> have been 
used in the transfer matrix computation: a half gaussian distribution of the angle tp (see text) and 
a flat one. 
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